Spin wave instabilities and field induced transitions in heavy fermions 
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We study phase transitions in heavy fermion systems due to spin-wave instabilities. One moti- 
vation is to determine the changes in the spin-wave parameters of a magnetically ordered heavy 
fermion system as it approaches a quantum critical point (QCP) by applying pressure. The other 
more actual is to provide an alternative approach, based on spin- wave instabilities, for the magnetic 
field induced transitions recently observed in antiferromagnetic heavy fermion materials. 
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INTRODUCTION 



Pressure and magnetic field have been used as exter- 
nal parameters to drive heavy fermion materials through 
different quantum critical points The theory of the 
heavy fermion quantum criticality has been developed 
starting on the paramagnetic side of their phase diagram 
0. In the ordered magnetic phase, usually with anti- 
ferromagnetic (AF) long range order, most studies have 
concentrated in determining the shape of the critical line 
as a function of pressure (P) and magnetic field (H) , 
in particular when Tn — ► |y. This gives information 
on the important shift exponent which is determined by 
the critical exponents associated with the QCP. However, 
in the broken symmetry phase, the magnetic excitations 
are spin-waves and a relevant question concerns what 
happens to these modes as the quantum critical point 
(QCP) is approached from the ordered side More 
specifically, we ask what happens to the spin-wave pa- 
rameters, its gap and stiffness as, for example, Ijy(P) 
goes to zero, i.e., P approaches Pc the critical pressure? 
The most direct way to answer this question would be 
to perform neutron scattering experiments, but these are 
non trivial. There are other less direct strategies to fol- 
low these modes as the system approaches a QCP. One of 
these is the electrical resistivity. Since magnons scatter 
the conduction electrons of the correlated metal, the tem- 
perature dependent resistivity has a component which is 
due to electron-magnon scattering. The temperature de- 
pendence of this component yields direct information on 
the spin- wave gap and stiffness Q. 

Recently, there have been many experimental studies 
on field induced transitions in magnetic heavy fermion 
systems Q . The analysis of these experiments takes into 
account correlated electrons with masses strongly renor- 
malized by the presence of a QCP at the critical field Hc2 
(see FigQ. However, field induced magnetic transitions 
are associated with soft spin-wave modes. Consequently, 
close to Hc2 these modes give an important contribution 
to the thermodynamic properties and electrical resistiv- 
ity, by scattering the quasi-particles. These contributions 



should be included in an analysis of this phase transition, 
besides that of other quasi-particles. 



,T x = (u/k B ) (H-H C2/ 



Paramagnetic 




FIG. 1: (Color online) Phase diagram of an anisotropic an- 
tiferromagnet in an external magnetic field in the easy axis 
direction^. The three phases, antiferromagnet, spin-flop and 
paramagnetic and the tricritical point (T) are shown. The full 
lines are lines of second order phase transitions. The dashed 
line is a line of first order transitions and the dotted line 
is a crossover line dividing the paramagnetic region in two 
regimes, as explained in the text. 



The magnetic field and temperature dependent phase 
diagram of an anisotropic antiferromagnet, with the ex- 
ternal field H applied in the easy axis direction is 
shown schematically in Fig^ General features of this 
phase diagram are the presence of a tricritical point and 
lines of first and second order transitions. Heavy fermions 
are in general strongly anisotropic and a relevant situa- 
tion for these materials is the metamagnetic case where 
there is no intermediate spin-flop phase at all and the 
system goes directly from the antiferromagnetic phase to 
the paramagnetic phase through a line of second order 
phase transitions, Hc2(T). 
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RESISTIVITY OF THE ANISOTROPIC 
ANTIFERROMAGNET IN ZERO MAGNETIC 
FIELD 

We start considering the case relevant for the study 
of heavy fermions systems, driven by pressure to a QCP. 
The spectrum of hydrodynamic spin-wave excitations in 
the antiferromagnetic phase, for H « Hci (T << T/v) 
(see FigOJ, is given by, 



u(k) = -fiH ± y/A 2 + Dk 2 

where the gap A is proportional to the anisotropy (cither 
exchange or single ion). In the isotropic case and in the 
absence of a magnetic field, the hydrodynamic magnons 
have a linear dispersion relation, to{k) = y/Dk, differently 
from those in a ferromagnet, for which, u>(k) = Dk 2 . 

We consider first a metallic antiferromagnet in a simple 
cubic lattice with no applied external field. Following 
Yamada and Takada 6] we write the electrical resistivity 
due to electron-magnon scattering as, 



p(T) = CI AF (H = 0, T) 



(1) 



where C is a constant. The quantity Iaf (0, T) is given 
by the following integral J6}, 



8 f°° k 2 



(2) 



where u(k) = VA 2 + Dk 2 is the dispersion relation of the 
antiferromagnetic magnons. In the limit B~ l = k B T « 
A, we deal with this integral changing variables and in- 
tegrating by parts to obtain, 



Iaf(0,T) 



2(k B T) 2 

£3/2 



ni 2 p-V0 



e yo(l-coshx) CQsh2 xdx 



yo(l — cosh x) 



where yo = A/k B T. The function F(x) — e 
such that, F(Q) = 1 and F(xf) = e _1 .The last condition 
implies yo(l — coshx/) = —1, or coshx/ = 1 + y^ 1 - This 
function cuts off the integral at x/ and we can rewrite 
the equation above as, 



Iaf(0,T) 
which yields, 



2(k B T) 2 

D 3/2 



vie 



Wo 



cosh 2xdx 



Iaf(P,T) 



A 2 

D 3/2 



e k B T sinh2a; 



(3) 



Since coshx/ = 1 + y 1 = 1 + k B T/A, sinh2x/ = 

- 1. Using that y^ 1 = 



2 ( 1 + % 1 )V( 1 + ^ 1 ) 

k B T/A « 1, we obtain 
2 V / 2A 2 



Iaf{0,T) 



D 3/2 



e k B T 



kuT 
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the first term in this expansion yields the result of Ref. |(| . 
Furthermore, it is easy to check that in the gapless case 
the temperature dependence of the resistivity due to lin- 
early dispersive magnons is, p{T) oc T 2 0. Eq. fallows 
to obtain from the measured p(T), the absolute value 
of the spin-wave gap. If p(T) is measured for different 
pressures it yields the trend of the stiffness D and gap 
A as a function of this parameter, for example, when 
approaching a QCP. In a recent analysis of a magnetic 
heavy fermion driven to a QCP by pressure, we obtained 
that A oc (Pc — P) , vanishing with pressure at the same 
critical pressure Pc that Tjv(P) vanishes while the stiff- 
ness D remains unrenormalized • It would be interest- 
ing to confirm these results by direct neutron scattering 
experiments 



ZERO TEMPERATURE FIELD INDUCED 
TRANSITION 

We now turn our attention to the next topic of this 
Report which deals with the phase transition induced in 
AF heavy fermions by an external magnetic field. In 
this case it is useful to start in the paramagnetic phase, 
with H > H C 2(T = 0) (see Fig. Q}. The results below 
apply in general even if the external field is applied per- 
pendicularly to the easy axis. The spin-wave spectra for 
H > Hq2 is ferromagnetic like, with a dispersion relation 
given by H@, 



;(q) = pH — S[J(0) - J(q)] 



(4) 



where the minus sign takes into account explicitly the an- 
tiferromagnetic nature of the exchange interactions (J). 
In a simple cubic lattice, the minimum of the spin-wave 
spectrum occurs for q = Q = (ir/a, 7r/a, it/a). For sim- 
plicity we have considered only the isotropic part of the 
exchange coupling 0, 0, 0|. Since ui(Q) = pH — 2SZJ, 
the condition lo{Q) = defines the critical field Hc2(T = 
0) = 2SZJ/p, below which the spin-wave energy be- 
comes negative signalling the entrance of the system in 
the AF or spin- flop phase. Above, Z is the number of 
nearest neighbors. Finally, writing q = Q + k and ex- 
panding for small k, the spin- wave dispersion relation can 
be obtained as, 



w(k) = p(H - H C2 ) + Dk 2 



(5) 



where the spin-wave stiffness at T = 0, D = SJa 2 (a is 
the lattice parameter). The temperature dependence of 
Hc'2 (T) is determined by that of the spin- wave stiffness 
D(T) This generally results in a critical line vary- 

ing with temperature as, H C 2(T) = H C 2(T = 0)-6T 3 / 2 . 

We can obtain the magnon contribution to the elec- 
trical resistivity in the presence of the external magnetic 
field as, p(T) oc If{FI, T), where in three dimensions (3d) 
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I F (H,T) = 



(ML 



C2 



3/2 



dy 



III) 



vVy - yo 

sinh 2 ^ 



(6) 



with y = A/k B T. The gap A = p(H - H C i) is the 
control parameter for the zero temperature field induced 
phase transition. 

We can distinguish two regions, separated by the 
crossover line T x = (p/k B )(H — Hci), in the phase 
diagram of Fig^ This line is fully determined by the 
value of the magnetic moment p, and by Hq2 ■ In region 
II, for T » T x and in particular for H = H C 2, the 
dominant contribution to the resistivity is due to gap- 
less magnons and given by, p cx (k B T / p,Hc2) 3 ^ 2 ■ In this 
region, for a fixed temperature T >> T x , the resistiv- 
ity as a function of magnetic field obtained from Eq[S] is 
shown in Fig'El ■ The resistivity has a cusp at Hc2- In 
fact dp/dH tx — (H — Hc2)~ 1 ^ 2 and diverges to — oo as 
H^H C 2 @- 

In region I, for k B T << p(H — Hc2) the resistivity is 
calculated as in the previous section and is given by, 



p cx 



k B T 
pH C 2 



3/2 



H, 



1 . 2 M (g 
3 k B T 



C2j 



e "b t (7) 
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FIG. 2: (Color online) Field dependence of the electron- 
magnon resistivity for a fixed temperature in region II of the 
phase diagram of Fig0 The inset shows the field dependence 
of the modulus of the derivative dp/dH, which diverges as a 
power law when H approaches Hc2- 

The thermodynamic properties of the system close to 
Hc2 can be easily obtained in the soft mode scenario. 
They are essentially those of a system of non-interacting 
bosons. The temperature and field dependent contribu- 
tion of the spin-wave excitations to the specific heat Cy 
can be shown to have a form similar to that of the elec- 
trical resistivity. We get, C v cx I F (H,T), with I F (H,T) 



given in Eq|B| Consequently, in region II, which includes 
the line H = Hc2,, the specific heat has a simple power 
law dependence, CV cx (k B T//iHc2) 3 ^ 2 - In region I, the 
field and temperature dependence of CV is the same as 
that of the resistivity given in Eq. [7| The magnetization 
at the critical field Hc2 decreases with a T 3 / 2 power law. 
In region I for H > Hq2 it decreases exponentially due 
to the gap, A = p(H - Hc%)- 



TWO DIMENSIONAL CASE 

Many heavy fermions systems due to their anisotropy 
in the crystalline structure have a regime close to 
their QCP where the critical fluctuations have a two- 
dimensional character Q . For this reason it is interesting 
to generalize the results obtained above for the case of 
two dimensions (2d). In 2d the relevant integral for the 
resistivity can be evaluated and we get, 



p2d(h,t) cx t 

where h = (H - 
limit (h/t) -> 0, 



2h , h 
— coth — 
t 2t 



4 In 2 - 4 In sinh ■ 



2t 



(8) 



Hc2)/Hc2 and t 
this reduces to 

0w(M)=/&(l 



k B T / pHc2- In the 



mi). 



(9) 



This resistivity varies nearly linear in temperature (Fig|3J) 
but for H — > Hc2 (h — > 0) it is not defined for finite tem- 
peratures. This is related to the fact that for an isotropic 
system there is no long range order in two dimensions at 
any finite temperature due to the creation of an infinite 
number of magnons. The systems we are considering are 
in fact three dimensional, although they may exhibit in 
some region of their phase diagram a 2d behavior. How- 
ever, this is just a regime which will cross over to 3d 
sufficiently close to the transition. Consequently this di- 
vergence of P2d(0,i) should not interfere with the actual 
observation of a nearly linear temperature dependence of 
the resistivity of real systems in some range of region II. 
In 2d the spin-wave contribution to the specific heat, as 
before, has a similar field and temperature dependence 
of the resistivity due to electron-magnon scattering. We 
find in this case ((h/t) — > 0) that in region II (Fig|2J), 

Cy(h,t) t 
cx a(l + ln~). 

The decrease of the magnetization with temperature in 
a fixed field in this region also has a field and temperature 
dependence similar to Eq. 

Finally, in region I the specific heat, magnetization and 
resistivity have an exponential term, as in Eq. Only 
the temperature dependence of the pre-factor is depen- 
dent on dimensionality. This makes it difficult to dis- 
tinguish a 2c? from a 3d regime in this region, either in 
thermodynamic or transport measurements. 
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FIG. 3: Semi-logarithmic plot of the resistivity p and specific 
heat Cv divided by T in 2d. Inset shows p and Cv versus T. 



CONCLUSIONS 



of the magnons and not the antiferromagnetic nature of 
the transition. In the extreme metamagnetic case of an 
Ising antiferromagnet in an external longitudinal field, a 
renormalization group approach yields a phase diagram 
where the T = transition at Hq2 is not associated with 
a fixed point. The point (H = Hc2,T = 0) is just an or- 
dinary point that flows to (H = 0, T = TV) as any point 
in the line of second order transitions H = Hc2 (T) , since 
the mag netic field turns out to be an irrelevant param- 
eter [llj. Notice that the spin- wave approach presented 
above is independent whether the T = transition at 
Hq2 is associated or not with a fixed point. 

We thank S. Bud'ko, J. Larrea and Magda Fontes 
for helpful discussions. Support from the Brazilian 
agencies CNPq and FAPERJ is gratefully acknowl- 
edged. Work partially supported by PRONEX/MCT 
and FAPERJ/Cicntista do Nosso Estado programs. 



We have obtained the electrical and thermodynamic 
properties of a metallic antiferromagnetic system due to 
spin- wave modes close to a zero temperature field induced 
transition. The existence of these excitations close to 
this phase transition is a general phenomenon, depending 
only on the appearance of a symmetry- broken state below 
Hc2- The resistivity is determined by the scattering of 
quasi-particles by ferromagnetic-like magnons which be- 
come gapless at the transition. The contribution of the 
spin-wave modes to the thermodynamic and transport 
properties is distinct in different regions of the phase dia- 
gram. In region I, this is exponentially suppressed by the 
gap A = fi(H — Hc2): but in region II these modes are ef- 
fectively gapless and give rise to important contributions. 
Notice that these contributions should be considered in 
addition to those of other quasi-particles. For example, 
in region I, it has been observed [4| a significant T 2 term 
in the resistivity probably due to electron-electron inter- 
action. This temperature dependence of the resistivity 
is different from that arising from spin-wave scattering 
which in turn is exponentially suppressed by the gap. In 
region II however, the spin-wave term is large and its 
power law temperature dependence is similar to that as- 
sociated with non-Fermi liquid behavior close to a QCP 
[Io| . This will make it difficult to separate this contribu- 
tion from that of other quasi-particles. The interaction 
of the latter with the spin- wave modes may even give rise 
to substantial mass renormalization. The same remarks 
apply for the thermodynamic properties in both regions. 

The temperature dependence of the thermodynamic 
properties at the critical field Hc2{T = 0) is similar 
to that predicted by the scaling theory of a QCP with 
a dynamic exponent z — 2 [lflj. In the present case, 
this value of z reflects the ferromagnetic-like dispersion 
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